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Divergence exponents of the first-order quantum correction of a two-dimensional hard-sphere Bose 
atoms are obtained by an effective field theory method. The first-order correction to the ground- 
state energy density with respect to the zeroth-order is given by £i/£q ~ \D — 2 "j hi7| a , where 
^q' D is the spatial dimension, and 7 is the gas parameter (7 = na D ). As D — > 2, a = a' = f . We show 

that the first-order quantum correction of the energy density is not perturbative in low dimensions 
of D < 2.2 regardless of any gas parameter which is much less that f . 

' PACS numbers: 05.30.Jp, 21.60.Fw, 03.70.+k 
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I. INTRODUCTION 



Weakly interacting uniform Bose gas is a fundamental topic of many-body Bose systems and has been studied for 
more than five decades. Bose-Einstein condensation(BEC) is known to be the only phase transition that does not 
^ require interaction. However, in the real system an interaction exists even in very dilute Bose gases, and often the 
diluteness is explained by a D-dimensional gas parameter: 7(1?) = naP <C 1, where n is the D-dimensional density 
and a is the s-wave scattering length. At zero temperature the chemical potential, energy density, number density, 
C ' and speed of sound are written as an expansion in powers of the gas parameter 7. 

, A fundamental approximation method of a weakly interacting Bose system was introduced by Bogoliubov in 1947 
I ■ In three-dimensions(3D) the ground-state energy density is given by 

o : 

O, £ + £ 1 + --- = £ [l + ciV7 + C27ln7 + ---]- (1) 



Similarly, the number density in 3D is given by 



no + m + ■ ■ ■ = no 



1 + v^+ 



(2) 
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£0 — 2irh 2 an 2 /m is the zeroth-order or mean-field energy density. The first-order quantum correction c\ — 128/15y / 7r 
was first obtained by Lee and Yang in 1957 for a hard-sphere Bose gas 0]. Later, it has been shown by Brueckner 
et al. [1, H[ that above results are generally true for any short-range potential with scattering length a. The second- 
order quantum correction ci = 8(47r — 3-\/3) /Sy^ that was originated from a three-body interaction was obtained 
by Wu et al. in 1959 [1, H, 0|- Recently, an effective field theory(EFT) method has been suggested to obtain the 
coefficients(Braaten and Nieto [1]). This method was introduced by Georgi et al. in 1990s [9I lid [Tlj. 

In two-dimensions(2D) an uniform Bose gas is governed by strong long-range fluctuations. The fluctuations inhibit 
the formation of a true long-range order. Therefore, 2D uniform Bose gas does not undergo BEC transition at finite 
temperatures [12j. However, this 2D system turns superfluid below a certain temperature while under confinement. 



The ground-state energy density of hard-sphere bosons in 2D has been studied by Schick in 1971 [1 31 ] . 

£0 + + — =£0 [1 + 0(1/ In 7)], (3) 

and number density is 

m + ni-i =7i 1 + : ,] , . H ■ (4) 

ln(l/ 7 ) 
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So = 2Trh 2 n 2 / (m\\nj\) is the mean-field energy density. The first-order number density in 2D is finite as ni/no = 
l/ln(l/7) ^ 1) but the first-order quantum correction of the ground-state energy density has not been obtained yet 
because an infrared divergence of an integral prevents a direct calculation. It has the divergence form of 



\D - 2\<*\lnj\°>" V ; 

where the exponents a and a' are positive constants. 

In this paper, we will calculate the two exponents of divergence a and a' analytically. The EFT has an advantage 
that it works even in non-integer dimensions. Therefore, we can rewrite the EFT in non-integer dimensions between 
two and three. Then, taking the limit of D = lim e _^ + 2(1 + e), we will obtain the two exponents of divergence a and 
a' . Furthermore, since the first-order quantum correction of the ground-state energy density converges in D=3 but 
diverges in D=2, we may find the dimensional border of the perturbative method which is effective between two and 
three. 

The paper is organized as follows: In Section II, we review the basic structure of EFT of a dilute Bose system in 
non-integer dimensions between two and three [H, Q EH • In Section III, we reproduce the zeroth-order or mean- field 
results in D-dimcnsions. In section VI, we obtain the first-order quantum correction in 2D and two exponents a and 
a'. Also we interpolate the first-order quantum depletion between 2D and 3D. Finally, in section V we summarize 
the results. 




II. EFFECTIVE FIELD THEORY BY AN EFFECTIVE LAGRANGIAN 



EFT is a general approach that can be used to analyze the low energy behavior of a physical system and takes 
advantage of the separation of scales to make model- independent predictions [tj Ell EH • The effective lagrangian that 
describes the low energy physics is written in terms of only the long- wavelength degrees of freedom that includes every 
non-renormalizable interactions, but to a certain order in a low energy expansion, only a finite number of operators 
contribute to a physical quantity ■ 

An uniform Bose gas can be described by field theoretic method with Hamiltonian TL — /iA, where (i is the 
chemical potential and A is the number operator. The number of atoms are conserved by the phase transformation: 
4>{Y,t) — > e tp(r,t). The free energy in ground-state of the system is 

F( l i) = {H- l iN) fl , (6) 

where (• • -) M denotes the expectation value in the ground state. Therefore, (7i) M = £(^) an d (A) M = n(fj,) are the 
energy and number in the ground state. The free energy T is given by all connected vacuum diagrams that are 
Feynman diagrams with no external legs |8j. The number density n is given by the expectation value of ('0''0) in the 
ground-state (we set h = 1 for convenience) 

n ( M ) = I J V^Vi>(^^)e lS[ ^^ ] (7) 

where and ip are complex field operators of bosons, S is the action from the Lagrangian Q, and Z is the grand- 
canonical partition function given by 

V^V^e lS{ ^ M . (8) 



The grand-canonical ensemble of uniform Bose particles with a two-body interaction is governed by the action 



S\tfA\ = jdtj d°X ^il + ^.+^-lg^^j 



(9) 



where m is the atomic mass, and g is the D-dimensional coupling constant which contains pairwise interaction between 
atoms. In D-dimensions g is given by [l6j 

4n D/2 a D ~ 2 1 

g(D) = — p: ^ , (10) 

yy ' m 2 2 - D r(i- f ) 7 °/ 2 - 1 +r(f -1) ^ ; 

where T is the Gamma function. It satisfies the two limiting cases. In 3D, it has the well-know form 



(11) 
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In 2D, it has the logarithmic form [171 1 181] 

When the temperature of a Bose system falls below the condensation temperature T c , we can write the quantum 
field ip in terms of a time-independent condensate <p and a quantum fluctuation field ip 

i> = (t> + i>. (13) 

The fluctuation field ip can be conveniently written in terms of two real fields ipi and ip 2 

$=±±j&. (14 ) 

In the uniform system <f> is the condensate order parameter and a real constant. It corresponds to a breaking of the 
global U(l) symmetry and takes into account the Bogoluibov shift of the field operator. 

The conservation of particle numbers requires that ip has non-zero momentum component so that (ip) = ef> and 
(ip) = 0. Therefore, the condensate order parameter <fi defines the density of condensed particles while ip defines the 
density of uncondensed particles. Then, the zeroth- and first-order quantum correction to the density are given by 

n = 4> 2 , m = ffiip). (15) 

Substituting Eqs. (fl~3f and (fT"4|) into Eq. ([9]) of the action S, we decompose the real part of the action into three 
parts BQEBl: 

S[<f>,1pl,1p 2 ] = Sdas[(f>} + Sf ref ,\<j>, 1pl,tp2] + Sint[<f>,ll>l,lh]- (16) 

S c i as is the classical part of the action. It does not contain any filed operator: 

S c i as [4>] = j dt j d°x {^4> 2 - i# 4 ^ . (17) 
Sf ree is the free part of the action. It is quadratic in ip\ and ip2- 



Sf rB e[<t>,i>l,lp2\ = I dt I d° x 



(18) 



where ip = dip/dt, and the two new variable X and Y are 



X = M-3 5 2 , (19) 
Y = fi-gcP 2 . (20) 



The terms 3g<fi 2 and g<p 2 in the X and Y are mean-field self-energies of the system. 
Sint is the interaction part of the action. It is the remaining terms: 



(21) 

J = <P(H- gtf) = 4>Y, (22) 



S in t[<j>,tpl,tp2] — j dt j d D X 

where 
and 

K = (23) 

The free part of the action in Eq. (|18p gives rise to a propagator, which can be used in perturbative framework. If 
we take the Fourier Transform in a momentum space, the fluctuating part is given by 

tp(t,r) = — ^2 5Z^(^n,p)exp{-iw n t + ip-r}, (24) 

V 71 — — OO p 

^ V - % . j M a. , . «• , Matsubara fr , w XW.e, t be „ < t be bee 



V 

action Sf ree is written as 



„, , 1 / £ Y —iui 



(25) 



x + xx + xxx + 
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FIG. 1: Loop diagram of the two-body interaction. The zeroth, first, and second loop from the left to right 



with the dispersion relation 




This is a general form of the dispersion and includes every information of the energy spectrum in the two self-terms 
X and Y. It is possible to diagonalize the propagator matrix in Eq. (|25l) by a field redefinition or renormalization, 
which is equivalent to the Bogoliubov transformation in the operator method. However, such a redefinition makes 
the interaction terms in the action more complicated and increases the number of diagrams. Therefore, we use the 
propagator matrix with off-diagonal elements to minimize the number of diagrams. 



III. THERMODYNAMIC POTENTIAL AND ZEROTH-ORDER RESULTS 



It is convenient to introduce the thermodynamic potential 0(/i, 0). The thermodynamic potential contains the 
information required to determine all of the thermodynamic functions. The free energy can be obtained by 

evaluating Q at a particular value of 4>. 

The sum of the vacuum graphs is independent of the arbitrary background condensate <p. Thus, the sum of connected 
vacuum diagram reduces to the sum of one-particle irreducible vacuum diagrams 

Fiji) = Clin, <M*)) (27) 
From Eq. © the ground-state energy density is written by 

S(ji) = J r (tx) + nfx. (28) 
By differentiating the free energy, we obtain the density and chemical potential, too. 

n(fJ,) = -—, n(n) = ±-. (29) 
a /j, an 

The n-loop contribution to the fl is denoted by fi n in FIG. 1. It is given by all one-particle irreducible vacuum 
diagrams and can be expanded in the number of loops: 

n(n, <t>) = n 4>) + Qi (ji, (f>) + n 2 Qu, 4>) h — . (30) 

If f2 is evaluated at a value of the condensate, all one-particle reducible diagrams vanish. Then, the free energy density 
in Eq. $T7$ is 

=n o (/Lt,0) + ni(Ai J <?i) + n 2 (M ! </>) + ••■ . (3i) 

Using Eqs. (fT5)l and (|14p. the condition of the condensate reduces to (ipi) = (1P2) = 0. This condition is equivalent 

to 

gn(/vji>) _ an o (^,0) ani(/x,0) an 2 (/j,^) _ , . 

d$ ~ d<j> + d<t> + d4> + • ' • - u M 

The number density in Eq. (|29p is obtained from the thermodynamic potential, too. 

»(M) = -f- (33) 

The free energy can also be expanded in powers of quantum corrections around the mean-field value 

T^) = To (At) + Tx (m) + & (fi) + ■ ■ ■ ■ (34) 

The loop expansion in Eq. (|31[) does not coincide with the above expansion of in powers of quantum corrections 
because of its independence of <j). To obtain the expansion of T in powers of quantum corrections, we must expand 
the condensate <f> around its classical minimum 0Oj which satisfies 

dn (fi,(t> ) _ 
84, ~ U ' 
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By expanding Eq. (|32[) in powers of <j) — </>o, and solving for 0, we obtain the quantum expansion for the condensate: 

<j) = 0o + 4>i + <f>2 + ■ ■ ■ , (36) 

where <p n is the nth-order quantum correction. For instance, the first-order quantum correction for the condensate 
4>i is obtained from the expanding of Eq. (|35p around <fio. The nonzero terms are 



dQoifJ-, <t>) d 2 {l (ii,<t>o) , . ar2i(^,0 o ) . „ 
_ ^ — 0oJ H — 1 = 0- 



d(f> 



Then, 04 in Eq. (|36[) is written as 



9fii(/x,$o) 



dm 



<9 2 fi (/x, O ) 



<90 2 



(37) 



(38) 



Therefore, keeping up to the terms of the second order and 
density is 



ii, the first three terms of the free energy 



Fo{n) = O (/x, O ), 

^(/X) = fii(/l,0 O ), 



<9fii(/z,0o) , , l9 2 rJ («,^o) 



prS'eias = -/X0 2 + ^30 4 - 



90 r 2 90 2 r1 ' 
The mean-field thermodynamic potential is given by the terms in the classical action in Eq. (|17[) as 

fl (fj,, <j>) 

From Eq. (|35p we obtain the classical minimum 0o as 

4>o = \ — 
V 9 

Therefore, the mean-field free energy density in D-dimensions is obtained from Eq. (|39[) 

,2 



.Fo(aO = ^o(/x, 0o ) 



V 



The mean-field number density is obtained from Eq. ((7|) as 

n (^o) = 77 — 1 = — 

dn g 

The mean-field ground-state energy density in D-dimensions is obtained from Eq. 



S, 



(D) _ A*o 



+ nfi = -g{D)n 



2g ' ^ u 2" 

In 3D, since g(3) = Aita/m from Eq. (jlip . the mean-field energy density is 

2nan 2 



(3) 



In 2D, since g (2) = 47r/mln(l/7) from Eq. JT2|) . the mean-field energy density is 

( 2 ) _ 2nan 2 
t o — 



mln(l/7) 

Therefore, this method reproduced the well-known mean-field energy densities in Eqs. ([T]) and successfully. 



(39) 
(40) 

(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 
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IV. FIRST-ORDER QUANTUM CORRECTION 

The first-order quantum correction of the ground-state energy is obtained in the following way. Substituting Eq. 
(|4"3")l into Eqs. (fTO|) and l|20p. we obtain the two variables X and Y at the minimum of <j>Q. 

X = n - 3g<ft = -2/i, (49) 
Y = fi - g<& = 0. (50) 

The propagator and the dispersion relation in Eqs. (|25[) and (|26p becomes 

(51) 

and from Eq. (|45|) 

p'2 / p2 



*=^\fa +2 *y (52) 

These are the original Bogoliubov results. E p is gapless and is linear for small wave-vectors. For large wave-vectors, 
the dispersion relation becomes E p ~ p 2 /2m + 2 / u ~ p 2 /2m + 2g{D)n. Therefore, the 2g(D)n represents the mean-field 
energy in D-dimensions due to interaction with the condensed particles. 

In the Bogoliubov approximation one makes a pair approximation to the Hamiltonian by neglecting terms with 
three and four operators On the other hand, in the Beliaev approximation, one goes one step further by calculating 
the leading quantum corrections to the quasi-particle spectrum [4|. This is done by including all one- loop diagrams. 

The free energy can be written with the inverse of the propagator. From Eq. © 



Z = exp|i j dt J d D xdetG 1 (w,p)j. 



(53) 



Therefore, 

T = TrlnG _1 (cj,p). (54) 

Note that IndetA = TrlnA for any matrix A. The one-loop contribution to the thermodynamic potential in D- 
dimensions, f^i = i In is obtained by the following way. 

dio f d D p 
2 J 2^ J (2ir) D 



ViM) = ~\ J ==; I T^lndetG^^)- 1 
1 f d D p 



2 J (2ir)°^ p 
From Eq. (|40p , we have the free energy component 



E p . (55) 



= SIi(jj.,4>o) = j—I ,-i(k), (56) 



where k = 4m/j = imgn and Iij are loop integral function in D-dimensions, which is defined as [1, [TJ 



>■' W J (2l) D p)(p 2 + K))/2' (57) 



where i and j are integers. This I^j is expressed by Gamma function, too 



r « (4n)"/* r(#)r(|) 



k 5 . (58) 



Note that in D=3. = k 5/2 /15tt 2 , 7^(«) = k 3/2 /3tt 2 , and /i 3 ^ (k) = -k 3/2 /6tt 2 . It has the useful relation 

of the first derivative 

^W = 4C(4 (59) 



Ground-state energy density is expressed as 

~ 2 ,y '" ' 4m 



So + Si + ■ ■ ■ = \gn 2 + i4?i(«) + ' ' ' • (6°) 
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FIG. 2: The first-order quantum correction of the ground-state energy density in D-dimensions in the logarithmic scale. The 
divergence in D—2 is clear. 7 itself is a function of the dimensions. 



Therefore, the contribution of the first-order quantum correction in D-dimensions is obtained as 



Si 
S 



r(-¥)r(¥! 

7r(o+i)/2r (f ) 



{mg(D)} 



D/2 D/2-1 



r (§) {2*-°t (1 - fj 7^/2-1 + r (f - i)} D/2 ' 



(61) 



Note that 7 = (no + rii H )a £> . 

Substituting D=3, we obtain the well-known first-order quantum correction of the ground-state energy density in 
3D as 

to 15v?r 

Taking the limit of D = lim e ^ + 2(1 + e), we obtain the first-order quantum correction of the ground-state energy 
density in 2D as 

£1 _ _ 7 £ 
£0 ~ 



> 1. 



-7 e + l e lii(l/7) 

For the singularities of the T function at non-positive argument, we used the following expansions 

(-iy 



(63) 



r(-n + e) 



+ h{n + l)+0(e) 



(64) 



where h(n + l) = 1 + 1/2 + 1/3 + - • -+l/n — 7^ and 7^ = 0.5772... is the Euler constant. In particular, T(e) = I/e—'Je- 
Also, we used the approximations: 7 e ~ 1 + elii7. 

In 2D there is a divergence, and we can obtain the two exponents of the divergence in Eq. J5| as a = a' = 1. We 
plotted Si/£q in FIG. 2 as the functions of the dimensions and gas parameter in the logarithmic scale. From the FIG. 
2 it is clear that the expansion in Eq. ([3]) satisfy the perturbative condition of £i/£o -C 1 only when D > 2.2. In low 
dimensions of D < 2.2, the first-order ground-state energy density is not perturbative regardless of the magnitude of 
the gas parameter. 

The number density is obtained from Eqs. (p~3j)-(fl~5ll as 



n = (f> 2 + ]r(tpi +4>l)- 



(65) 



Therefore, taking the one-loop effects into account 

dui dPp 

no + Hi ' 



p 2 /2m + 2mE 2 /p 2 
co 2 -E 2 + iE p 



U0+ 2j 2n(2*)° 



n + ~ h t \{4mgn) + ^I-i-i{Amgn) 



n 



! + r (-f ) l 2r (zS zEtzr)} { )D/ 2(d-2)/2 

2 3 7r (D+2)/2 r (D\ V yj 



(66) 



8 



0.12 

o 




1 ' ' ' ' 1 

3 2.8 2.6 2.4 2.2 2 

Dimension, D 



FIG. 3: The first-order quantum correction of the number density at 7 = 10 4 as a function the dimensions. 



Then, the first-order quantum correction of the number density in D-dimensions is obtained as 

21 = ^ 2 - 2 "- 2)/4 r (-f) {2r (s±i) -ra 

»o 23 ~ D r(f){22-^r(i-f) 7 W2-i + r (f -i)} D/2 

In 3D we obtain the well-known first-order quantum depletion as 

^ = ^77- (68) 
no 3vi" 

In 2D taking the limit of D = lim e ^ + 2(1 + e), we obtain the first-order quantum depletion as 

ni ^ 7 £ 

1 



ln(l/7)' 



(69) 



We used the approximations: T(— -2) ~ 1/e and 2r( i ^±I) — r(^i) ~ 2y / 7T£ . This result is the same as the Schick's 
value in Eq. We plotted ni/no as a function of dimensions and gas parameter in Fig. 3. It is increasing 

monotonically as the dimension is decreasing from 3 to 2. The depletion is stronger in low dimensions. 



V. SUMMARY 



EFT is a theory of symmetries. Once the symmetries have been identified, one writes down the most general local 
effective Lagrangian consistent with these symmetries. At zero temperature, the symmetries are Galilean invariance, 
time-reversal symmetries, and the global phase symmetry. These symmetries restricts the possible terms in the 
effective action. 

The first-order quantum correction to the ground-state energy density of a 2D uniform Bose atoms has the form 
of the divergence: £i/£o ~ \D — 2| _Q | ln7|~ a . Applying an effective field theory method to the hard-sphere boson 
system in D-dimensions, we obtained the D-dimensional free energy and thermodynamic potential up to the one-loop 
result. This general form in D-dimensions reproduced every well-known result in 3D and 2D. 

Then, taking the limit of D = lim e _ +0 + 2(1 + e), we obtained the zeroth- and first-order quantum correction in 2D. 
Finally, we obtained the two exponents of divergence a and a' analytically as a = a' = 1. We also showed that the 
perturbative expression of the ground-state energy density is effective only in low dimensions of D < 2.2. 
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